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This paper describes the effects obtained in phase contrast imagery with partially coherent illumination 
of phase objects. The objects chosen for study include simple Zernike type of periodic phase gratings of the 

form (i) A(u) = l-\-i- A-cos mi and (ii) A(u) = l-{-t* ^J A n -cos nom, with axial or oblique partially coherent 

71=1 

illumination. The concept of "effective source" has been utilized, and the nonlinearity in the intensity 
transmission due to partial coherence has been analyzed in terms of "fundamental contrast" and "remnant." 

1. Introduction 

1 1 is known that ordinary methods of observation are not suitable for obtaining information 
about phase objects. Special techniques have to be used to render phase variations visible, 
namely, the dark-ground method of observation in which the direct beam illuminating the 
object is stopped. There is the other method, the Schlieren method, which differs from the 
above in that not only the direct light but also a part of the diffracted beam is also stopped 
from entering the ultimate detector. However, the most powerful method in this respect is 
due to Zernike [1934, 1942] which makes use of a phase modulation of the central or zero order 
of the diffraction produced by the object. In fact, it was shown by Zernike that with this 
method, which he termed as "phase-contrast" method, phase variations in the object are 
practically linearly transformed into changes in intensity and thus made observable. Later, 
Hopkins [1952] pointed out that because of the finite size of the aperture, some details of phase 
structures could be seen in ordinary microscopes as well. 

The Zernike theory of the phase contrast method basically makes use of the Abbe theory 
of microscopic imagery in coherent light. We shall then refer to this as the Abbe-Zernike 
theory. A simple and obvious extension may be made to the theory to include the effect of 
amplitude modulation of the diffraction spectra, thus including the dark-ground and Schlieren 
methods of observation as well. In the original elementary theory of Zernike there were 
certain implicit assumptions which were later treated by Picht [1936] and Kahn [1955]. Oster- 
berg [1944] has also made a thorough study of the phase-contrast method as a problem in 
diffraction, in which he showed that the phase-plate characteristics could simply be treated 
as a "coating function" modifying the pupil function of the microscope objective. His basic 
diffraction integral, which he terms as "generalized transport integral" includes the effect of 
not fully coherent illumination, though in the final analysis he restricts himself to a narrow 
cone, or effectively coherent illumination of the object. 

In the meantime, the concept of partial coherence has been utilized by Hopkins [1953] 
in formulating what he describes as the generalized Abbe theory, the generalization being due 
to his inclusion of the effect of partially coherent illumination of the object. Steel [1958] 
also gave a generalized treatment of image formation in partially coherent light. We followed 
up the obvious next step in generalizing the Abbe-Zernike theory of "phase and amplitude 
contrast microscopy in partially coherent light." As a first step in this direction De and 
Som [1963] extended the study of Hopkins [1952] to show that phase structures could be ob- 
served in ordinary microscopes even in partially coherent light, the extent of observation 
depending on the magnitude of phase variations and the degree of partial coherence. 
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In this paper, we now discuss the effect of phase and amplitude contrast technique on 
the images of simple Zernike-type objects represented by A(u) = l J r i A(u), in which the mag- 
nitude of A is small enough to justify the presumption that the object shows no variable 
intensity transmission. In our treatment, we have utilized the theorem of effective source 
as postulated by Hopkins [1953, 1957], and the observation of Osterberg that the phase and 
amplitude contrast plate characteristics may be simply associated with the pupil function of 
the objective. This treatment may be further extended to cover a general type of phase 
grating represented by A(u)=exip [i A(u)], of which the Zernike-type object is just a special 
case. We may further assert that our method is perfectly amenable to the consideration of 
even nonperiodic objects in a rather straightforward manner, and also of the effects of small 
residual aberrations in the objectives. Slansky [1962] also studied a similar problem, though 
his treatment is restricted to the case of Zernike-type objects only. 

2. General Solution 

In deriving the general solution of the problems of diffraction in phase and amplitude 
contrast microscopy in partially coherent light, we shall make use of the formulation by Hop- 
kins [1953]. In his so-called generalized Abbe theory in partially coherent light, Hopkins 
showed that, in all practical cases, the illumination of the object may be considered as due to 
a self-luminous source placed at the condenser exit pupil. This is termed the "effective source/' 
and insofar as the degree of partial coherence on the object plane is concerned, there seems to 
be no difference of practical significance between the two types of illuminating systems — • 
viz, the Kohler and the critical illumination. This and the further consideration of the 
progress of mutual intensity between a pair of points through the object to the objective aper- 
ture shows that this "effective source" may be orthoscopically projected on to the objective 
aperture through the object point. 

Earlier Osterberg and Wilkins [1949] showed that with objectives satisfying Abbe's sine 
condition, the diffraction plate characteristics may be simply treated as modifying the pupil 
function of the objective. Thus, the overall pupil function of the phase- and amplitude-con- 
trast objective may be taken as 

J(x,y)=j(x, y)p(x,y) 

where p(x, y) is the complex amplitude transmission of the diffraction plate. That this is pos- 
sible is seen from the fact that (x, y) denotes not only the coordinates of the point on the pupil but 
also the angles a ray makes with the axis as it comes from the object point on the axis to the 
point on the pupil. Thus, this also refers to a corresponding point on the back focal plane on 
which the diffraction plate is situated. 

Based on these two premises, the final intensity distribution in the image plane is given by 

B'{u\ v')= j I y(x , y )\4>(x 0j y ; u', v')\ 2 dx dy (1) 



where 



u =----n -sin a •£ , 

A 



2tt 

z X 



n • sin a • r\ , 



Xo ~h' y °~h' 

and where 

u f v' — reduced coordinates of a point on the image plane, 
£', 7)' — coordinates of a point on the image plane, 
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#o> Vo — fractional coordinates of a point on the pupil, 

sin a — semiangular aperture of the objective referred to the object point, 
&o? b — coordinates of a point on the pupil or on the effective source, 
h — radius of the pupil, 

n' — refractive index of the medium between the object and the objective, 
y(%o, 2/o) — i« the effective source function described in terms of the pupil coordinates, 
\<t>(x , y ; u', v')\ 2 — partial intensity at the point (u', v') of the image plane due to an element 
dxtflyo of the effective source. 

The physical significance of the above integral is obvious. To obtain the partial intensity, 
\<t>\ 2 , we may proceed as follows: 

Let A(u, v)= complex amplitude transmission of the object, which may be Fourier analyzed 
into 

a(Xy y) = \\ A ( u > v)e- iiux+vy) dudv (2) 

giving the components of the spectrum of A(u, v). 

If we consider the object to be illuminated by an element dx () dy of the effective source at 
x , y , the spectrum would be shifted by a corresponding amount. Thus, the amplitude 
distribution of the spectrum would be given by a(x—x , y—yo), which, after passing through 
the objective and the diffraction plate, will be modified to a(x— x , y—yo)j(x, y). The partial 
intensity, then, is given by the squared modulus of the inverse transform of the above. In 
other words, 



*< 



*6, Vo) W, v')= [ja(x-x 0} y-]h)f (x, y)e Hu ' x + v ' v) dxdy (3) 



gives the complex amplitude distribution in the image plane due to an element of the source 
at x , yo. It may also be interpreted as the complex amplitude distribution in the image plane 
when the object is illuminated by a plane wave front proceeding in the direction making an 
angle x , y with the axis. 

It now remains to evaluate explicitly the function pertaining to a specific object, which 
is possible if the composite pupil function is known. 

3. Solutions of Cases of Periodic Grating Objects of the Zernike Type 

We are now in a position to obtain solutions of specific cases. First, we shall consider 
those types of diffraction (or Zernike) plates which have a phase and amplitude contrast 
coating over an area which has the same relative size as the effective source. The complex 
amplitude transmission of this area, which we shall call the "conjugate area," is, say, given 
by he** j while that of the remaining area, the so-called "complementary area," is unity, h in 
the above implies amplitude contrast and \f/, phase contrast. 

Let us first consider the simple Zernike-type object: 

A(u) = l+i>A(u) (4) 

in which A is so small that A 2 is neglected. In other words, the intensity transmission of the 
above object is constant while the phase varies as A(u). Assuming A(w) to be periodic, we 
have 

A(^)=S A n • exp [ — inoou] (5) 

where 



1 f+% 
A n =— A(u) exp [inam]du 

* tfo J ~ u () 



2uq is the length of the fundamental period. co= (71-/1/ )- 
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The Fourier-spectrum of the object function (4) is thus given by 



+ o 



a(x, y)=d(x, y)+i S A n • 5(x+nu, y) (6) 

where 5(x, y) is the Dirac delta function. 

Thus, we have for <£ as defined in eq (3) 

(k(x ,y ',u f y)= I 6(x—x , y—y )+i S A n -8(x—Xo+na), y—y ) I/Or, y) exp {i(u'x+v'y) }dxdy 

=lf(^yo)+i^2^n-J(xQ—ncc,yQ) exp {—inuu'}] exp {i(u'x +v'y )} 
or 

l*l 2 =l/(»o, 2/o)| 2 +2-Im {/U, !/o)Sa:/*(x -^, y ) exp (iw ; ) } (7) 

in which the asterisk denotes the complex conjugate, and Im { } indicates that the 

imaginary part is to be taken. Terms involving A 2 have been neglected. The final intensity 
would then be obtained by an integration of |</>| 2 over the effective source y(x 0j y Q ) as indicated 
in eq (1) above. The integration involved is of the form 

O(0, n) = J j t (xo, yo)f(xo, y )f*(x —nw, y Q ) dx dy (8) 

which is the generalized partially coherent response as described by Hopkins. 

Since we shall be concerned with pupils or Zernike plates which are symmetrical about 
the y -2ixis, Jfao+uu, yo)=J(x Q —ucx), y ) and if, further, we restrict ourselves to cases where A n 
is real, then eq (7) reduces to the following simple form: 

°° _ _ 

M 2 = l/(z<>, Vo) | 2 +2 2 Im {j\x , y )f*(xo—nw, y ) (2A n ) cos na>u' . (7a) 

u = l 

The final intensity is thus given by 

B'(u')=l+2S (2A W ) Im ^ (0 :! ?) } cosnw' (9) 



n=l 



C(0, 0) 



in which a normalization is effected by the constant term (7(0, 0) . The (7's in the above equation 
are as defined in eq (8). The summation has been shown to extend from 1 to N, where 
Nu< 1+p since (7(0, n) shall be for larger values of n. 

The problem now reduces simply to an evaluation of the (7's with given Zernike plate 
and effective source characteristics. 

4. Evaluation of C(0,n) 

In this section we shall describe the mode of evaluation of (7(0, n). We shall consider 
first axial illumination of the object with a circular disk-shaped effective source. The radius 
of this source projected on the aperture of the objective is, say, p. The actual radius of the 
condenser aperture may be obtained from the following [viz: Som, 1963]: 

n c - sin a c =n • sin ao=n • p • sin a 

where n c and n are the refractive indices of the media between the condenser and the object, 
and between the object and the objective respectively, a c is the semiangular aperture of the 
condenser from the object point, a is the corresponding semiangular aperture of the effective 
source projected on the objective, and a is the semiangular aperture of the objective. In what 
follows we shall consider only cases in which p<0- 
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Second we shall consider annular illumination, i.e., cases of an annular condenser aperture 
with radii 1 and p a <0. The objective shall be assumed to be free from aberrations. Further, 
the conjugate area of the Zernike plate shall have the same size and shape as that of the "ef- 
fective source." 

(i) Axial illumination: 

Let us recall the expression for (7(0, n) in eq (8) 

C(Q,n)= I I y(x , y Q )J(x 0) yo)f*(x —no) ) yo)dx dy . 

As indicated above, we may write 

y(x , y ) = l within a circle 

= outside. 
fao, Vo)=J(x , y ) p(x ,y () ) 
j'(x , y ) = 1 within a circle 

= outside 
p(x 0} y )=he ixl/ within a circle 

= 1 outside. (10) 

Figure 1 shows the geometrical implication of the integral (8). Evidently the region of inte- 
gration is complicated by the various limits imposed on the functions that appear in the inte- 
grand. The values of (7(0, n) have been obtained by subdividing the area of integration into 
different areas of overlap of two circles of same or different radii. These latter areas could be 
easily determined analytically, and with the values of the function within these different 
subdivided areas, one would obtain the final values of (7(0, n) in terms of h, \p, n, co, p. They are 
given as follows: 

Im { C(0,n) }= / 1-D (— ) \ h sin $ ■ tt P 2 for co<l-p 

= | l-Z?o (y)-» Do ( ? f)+27 A(2 -^+7 2 ) y h sin + . 



7rp 2 



for 1 — p<w< Vl— / 



where 



{ ~ D ° (!f) + l D ° (^?) + 2? d ° (2 ■ no)+X2) } h sin * 



1-p 2 -iV , (l-p 2 )+A 2 f . , 2 



• 7Tp" 



C(0,0)=7T P 2 -A 2 (11) 

and 

A(aO=; [2 • cos" 1 ^-sin ■ (2 cos" 1 0]- (12) 

It should be noted that D (x) gives the response of a perfect lens system imaging a self- 
luminous object. These values are already tabulated, and hence the formulas for (7(0, n) given 
as above permit a straightforward evaluation, and also may be analyzed for aberrational effects 
with relative convenience. 
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Figure 1. Evaluation of C(0,n); axial illumination. 




Figure 2. Evaluation of C(0,n) ; annular illumination. 



(ii) Annular illumination: In this case, we have 

tC^o* 2/0) = 1 within an annulus bounded by circles xl+yl=l and #o+2/o=Pa 

=0 outside 
f&oy Vo) =1 within a circle xl-{-yl=l 

=0 outside 
p(%o> yo)=he i * within an annulus bounded by circles xl-\-y^=l and xl^-yl=p 2 a 

= 1 outside. 

Figure 2 shows a typical geometrical representation of the area of integration in this case. The 
values of C(Q, n) are as follows: 

lm{C(0,n)}=Trp 2 a \ 1— Z? ( — )Usin^ fornco<l— Pa 

L Vpa/ 2 \ nu Pa J 2p 2 a \ 



71 ca 



— ) \h sin \p 



for 1 — p a <ttW<A 1— Pa 



L \Pa/ 2 \ 71cop a / 2p a \ 7lCO /J 



aO,0) = 7r(l- Pa )A 2 . 

5. Some Results 



for n co > a^— Pa- 
(13) 



(i) Axial illumination: 

(a) Singly periodic grating: 



We shall first consider the results for a singly periodic phase-grating. Such a grating is 
represented by 

A(u) = l J ri -A • cos &u 
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Hence, from (9), we have 



/>"((//) = l-\-2ai ■ A ■ cos cot 



where 



«i 



Im [C(0, l)] 
<?(0, 0) 



(14) 



(15) 



With coherent illumination, i.e., in the limit of the effective source y(x Q} y ) tends towards 
the Dirac delta function at the origin, a x tends towards (sin xp/h) for all values of co<l, so that 
the intensity distribution is given by 



B' c (u') = l+2 • ^ • A • cos im\ 



(16) 



Then the contrast in the image with partially coherent phase and amplitude contrast relative 



to that with a limitingly narrow-cone axial illumination is given by a= 



sinxp 



Further, with 



narrow-cone axial illumination, the contrast will be zero for co>l, while in the partially coher- 
ent case, this limit is co>l + p when a becomes zero. Thus, there is at least a possibility that 
there would be a gain in resolving power of the phase-contrast objective with a finite width of 
the illuminating beam. However, this is obtainable at the cost of contrast. 

The values of a have been calculated for various values of co and p. Figure 3 shows the 
results graphically. The following summarizes some of the pertinent observations in respect 
of the above. 

(1) So long as 2p <co < 1 — p, the imagery in partially coherent light is indistinguishable from 
that with narrow-cone illumination in so far as the contrast is concerned. Hence, if one can 
have a prior knowledge of the range of frequencies in the object, it might be possible to choose 



J-OOr 




773-721—05- 



Figure 3. Contrast in the images of phase gratings; axial illumination. 
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a finite size of p to obtain a better photometric efficiency of the setup. Figure 4a shows the 
formation of the spectra under these conditions. All the spectra are wholly within the pupil 
and completely separated. That is why the ultimate image shows no difference in contrast 
from what is obtainable with coherent illumination. 

(2) If oo<2p < 1 — p, the spectra are still wholly within the pupil, but they mutually overlap 
to a certain extent. Thus, not only the entire zero-order spectrum, but also parts of the 
first-order spectra are affected by the phase plate. This latter neutralizes the phase modula- 
tion, and hence contributes to a general loss in contrast. Figure 4b represents the state of 
affairs under these conditions. 

(3) If l+p>co>l— p, parts of the first-order spectra fall outside the pupil and hence do 
not contribute to the image formation. The result is a further decrease in contrast. Figure 4c 
represents these conditions. 




Figure 4a. Formation of the spectra: 2p:<co:<l — p. 
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Figure 4b. Formation of the spectra; co<2p<l — p. 



FrGURE 4c. Formation of the spectra; l + p> «> 1 — p. 
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Figure 5. Useful frequency-range for partially 

coherent phase-contrast images. 




The values of a as shown in figure 3 may now be analyzed in the light of the above observa- 
tions. Figure 5 shows (i) the range of frequencies which shall have a contrast same as in the 
coherent case for different values of p, and (ii) the range of frequencies for which the relative 
contrast is equal to 0.8. This gives one an idea of the permissible maximum size of source 
that can be used for a satisfactory phase contrast image. 

The curves in figures 3 and 5 have been drawn for cases when h=\, that is, when no 
amplitude modulation is used. Further sin \p has been taken as unity as well. This implies 
that the phase of the zero-order spectrum is advanced or retarded by ir/2. Noting that a x is 
proportional to sin x//, the maximum phase contrast efficiency is obtained at this value. How- 
ever, «! is seen to be inversely proportional to h, meaning thereby that the contrast may be 
enhanced by using a slight amplitude contrast as well. This, of course, is possible at the cost 
of the mean level of illumination of the final image. 

(b) Square-topped phase-grating: 

This is the case of a simple periodic structure where troughs and crests represent differences 
in optical path only. 

Let 

A(u) = l+iA(u) (17) 

where A(u) is periodic in u with a fundamental length of period 2u . Let, over one such 
fundamental interval, 

A(u) = aito<\u\<u 



=A 



\u\ <au a . 



The Fourier spectrum is given by 



sin arnr . 
A n =A- n =^ A- 



(18) 



a in the above gives the ratio between the line width to line separations. Hence, from eq (9), 
we have 



where JVw <l+p. 



2?'0/) = l + 2A.f:2. S ^^^.r^g^l 
t=i nw L C(0, 0) J 



cos nuu 



(19) 
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In the limit of a point source, the expression for intensity will be 
B (, ) = 1 + 2A g (2 • ——y^-) cos 



nun 



(20) 



where iVco<l since Im [C(0, n)] in the coherent limit will be h sin \p for all values of nco<l. 

Figures 6a to d show the image of such gratings (A=0 • 5) for four different frequencies, and 
for various effective source sizes at each frequency. The dotted curves represent an ideal 
case of linear unitary transformation of phase variation into amplitude variation. 

The expression for B' (V) in eq (19) or for B' c (u') in eq (20) shows that the final image is 
given by a truncated Fourier series, and hence one would expect to observe the so-called Gibbs' 
phenomenon of large oscillations in the region of the edges of the square-topped grating. In 
the presence of harmonics in the image intensity distributions it is indeed a bit difficult to 
define a unique contrast for such cases. The presence of these harmonics is in fact an indica- 
tion of the effect of nonlinearity in the intensity transmission in partially coherent light, and 
as yet no unique specification has been suggested for a proper measure of such nonlinear effects. 
One suggestion is to express the contrast in terms of the fundamental frequency distribution, 
and lump together the effects of higher harmonics in a so-called "remnant" [Ingelstam, 1964]. 
This "remnant" may be described as the square root of the sum of the squares of all coefficients 
other than the constant and that of the fundamental frequency component. Figure 7 gives a 
plot of this fundamental contrast and the remnant against source sizes for different frequencies. 

It is of interest to note that the remnant always decreases with increasing frequency, 
falling gradually down to zero. In fact, if tV<2, there are no harmonics in the image intensity 
distribution. Thus, for a frequency up to 0.5, harmonic distortion would be present for all 
source sizes, though its magnitude as reflected in the value of the "remnant" may be negligibly 
small. For values of co lying between 0.5 and 1, the ''remnant" would be zero only if p is less 
than a certain value given by p<2co — 1. 




Figure 6a. Images of square-topped phase gratir* 
axial illumination; u = 0.25. 




-u)u / 



Figure 6b. Images of square-topped phase gratings; 
axial illumination; 03= 0.50. 
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Figure 6c. Image* of square-topped phase gratings; 
axial illumination; 03= 1.0. 



Figure 7. "Fundamental contrast 77 and 
u remnants ,J in the images of square- topped 
phase gratings; axial illumination. 




Figure 6d. Images of square-topped phase] \j rati ngs; 
axial illumination; 03= 1.5. 



hOr 




In the cases for which computation has been done, a has been chosen as equal to 1/2, im- 
plying equal widths for the troughs and crests of the phase grating. In these cases, the spectra 
consist of only the odd harmonic components and hence the "remnant" would be zero if A^<3; 
i.e., if p<3co— 1, co>0.33. This expresses the range of validity of using the criterion of "funda- 
mental contrast" for the study of the efficiency of phase-contrast devices in partially coherent 
light, 

(ii) Annular illumination: 

(a) Singly periodic grating. 

The expression for intensity is the same as in the previous case, except that the values of 
1m [C(0, n)\ to be used are as given in (13). 
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If the width of the annular source tends towards an infinitesimally small value, the area 
integral for (7(0, n) will tend towards a line integral, and we shall have 



lim Im[(7(0,7i)]=2.cos- 1 .('^AsiniA 



Pa^ 1 

lim [C(0,0)]=27rh 2 . 

Pa^ 1 



n) 



Thus, unlike the case of axial illumination discussed above, we have in this case a contrast 
varying with frequency even with a source of infinitesimally narrow width; we have, in this 
limiting case, 



cos" 



B' C =1+2A 



'(t*) 



COS coll 



(22) 



The difference with the corresponding expression (16) for the axial case is to be noted with 
interest. 

Evidently in this case the limiting resolvable frequency is no)<2. This is just a theoretical 
limit, as the contrast also tends towards zero at the limit. 

If we now consider the case of the annular source with a finite width, we shall find that 
Im [C(0, n)] becomes zero for 7ico>l + p fl which implies that we lose on the count of ultimate 
resolving power. Figure 8 shows graphically the values of a n =Im [C(0, n)]/C(0, 0) as a function 
of frequency co for various values of p a . The dashed curve in figure 8 indicates the range of 
frequencies which are so imaged that the relative contrast is nowhere less than 0.8. 

(b) Square-topped grating: 

As in the case of axial illumination, we have, in the limit of a line source coinciding with 
the margin of the pupil, 



where 



B c (u')=l+2A g (—^—) {-) (-f) cos W 



3 n =cos 



nw 



andN<(l-\- Pa )/w. 



(23) 




Figure 8. Contrast in the images of phase gratings; 
annular illumination. 
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With a finite width of the source, we have an expression similar to (20) above. Figures 9a to d 
show the images of such gratings (A=0.5) of a number of frequencies <o=— > where 2u is the 
length of the period, and a=| implies equal widths for the troughs and crests of path-differences 




Figure 9a. Images of square-topped phase gratings; 
annular illumination: w = 0.25. 




Figure 9c. Images of square-topped phase gratings; 
annular illumination; u=1.0. 
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Figure 9b. Images of square-topped phase gratings; 
annular illumination; 03 = 0.50. 




Figure 9d. Images of square-topped phase grating 
annular illumination; o>=1.5. 
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Figure 10. "Fundamental contrast" and "remnants" 
in the images of square-topped phase gratings; 
annular illumination. 



in the grating. The dotted curves represent an ideal case of linear unitary transformation of 
phase variation into amplitude variation. Here also we may make the same remarks as done 
in the case of axial illumination of square-topped gratings. Figure 10 describes the situation, 
and the conclusions are also similar. Remembering that the coherent limit of resolution in 
this case is 2, we find that practically the upper half of the range of resolvable frequencies is 
still assessable in terms of the "fundamental contrast" only. 



This study has been undertaken under the sponsorship of the National Bureau of Stand- 
ards. The authors are deeply indebted to the stimulating discussion with S. C. Som on the 
various aspects of this study. 
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